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We show theoretically that multipartite entanglement is generated on a massive scale in the spec- 
trum, or optical frequency comb, of a single optical parametric oscillator (OPO) emitting well above 
threshold. In this system, the quantum dynamics of the strongly depleted pump field are respon- 
sible for the onset of the entanglement by correlating the two-mode squeezed, bipartite-entangled 
pairs of OPO signal fields. (Such pairs are independent of one another in the undepleted, classical 
pump approximation.) We verify the multipartite nature of the entanglement by evaluating the 
van Loock-Furusawa criterion for a particular set of entanglement witnesses deduced from physical 
considerations. 



I. INTRODUCTION 

The generation of massively entangled states is of great importance for quantum information. For quantum commu- 
nication, good examples are multiparty quantum teleportation [1] and quantum secret sharing [2]. For measurement- 
based quantum computing, cluster states 4] are known to enable one-way quantum computing [l, (>]. Constructing 
large-scale quantum registers and processors is therefore one of the prime objectives of experimental quantum infor- 
mation, along with the suppression or alleviation of decoherence. 

In most cases, the approach to scaling up the size of quantum registers or processors is a "bottom-up" one, in 
which individual Qbits (following Mermin's more harmonious spelling [ ]) are put together to form, say, an entangled 
quantum register [8]. Now, there are, indeed, extremely few examples of "top-down" approaches to multipartite 
entanglement, in which a single physical system enables intrinsic generation of multipartite entanglement over a large 
scale. To the best of our knowledge, there are but two such systems. The first one is the individually trapped atoms in 
an optical lattice initially loaded with a Bose-Einstein condensate subsequently undergoing a Mott insulator transition 
[9]. The second one is the ensemble of entangled quantum modes of light, a.k.a. "Qmodes," defined by the resonant 
frequencies — or quantum optical frequency comb (QOFC) — of an optical parametric oscillator (OPO), in which the 
QOFC is entangled by the OPO's nonlinear crystal [10, 11]. Recently, the simultaneous generation of 15 identical 
quadripartite "square" cluster states was realized experimentally over 60 Qmodes of a single OPO \. 

The QOFC entanglement experiments mentioned above necessitate an exquisitely sophisticated OPO [12], operated 
below threshold [13], and in which two or three different nonlinear interactions must be phasematched [14, lO]. 

In this paper, we present the theoretical discovery of massive multipartite entanglement generation in a much 
simpler system and in a completely different regime. The system is but a standard OPO, in which only one nonlinear 
interaction is phasematched. In addition, the OPO must be operated well above threshold. It is somewhat surprising 
that such a simple, well-known system might lend itself to the generation of such an exotic quantum state as a 
massively multipartite one. In particular, we emphasize that the entangling interaction is only pairwise. It is the 
fact that all entangled pairs are derived from the same, strongly depleted pump field that generates the multipartite 
entanglement by way of a bona fide 3-field Hamiltonian. This is therefore a fundamentally different situation from 
that of the below-threshold OPO in which pairwise interactions are chained and all their pump fields are undepleted, 
yielding quadratic nonlinear interactions [id] in lieu of cubic ones. 

This paper is organized as follows. In Section 2, we introduce the system Hamiltonian, distinguishing between the 
depleted and undepleted pump cases. In Section 3, we solve the equations of motion for the system by employing a 
linearization procedure. We are certainly aware that more sophisticated treatments exist [13, 1(1-10] and may indeed 
be interesting to use in order to explore this system further. In particular, it is worth mentioning the new physics of 
noncritical squeezing generation — that is, squeezing independent of the system parameters such as pump amplitude — 
in the transverse spatial modes of an OPO, which was recently predicted via the phenomena of spontaneous symmetry 
breaking [20, 21] and pump clamping [22], the latter having already been observed in the laboratory [23]. Also, in 
this novel regime, the well-known, laser-like (and usually slow) phase diffusion process of an OPO [24, 25] becomes 
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entwined with the squeezed variables and affects detection [21], which isn't usually the case in a critically squeezing 
two-mode OPO [24, 25]. As the present paper doesn't pertain to noncritical squeezing, we have set aside this issue 
of phase diffusion for further studies, under the hypothesis that its effect may be similar to that in the usual critical 
squeezing situation. We therefore focus here on the nontrivial new results obtained from the simple approach adopted 
here. In Section 4, we use the multipartite inseparability criterion derived by van Loock and Furusawa [ : ] to establish 
the existence of multipartite entanglement in the optical frequency comb of a single OPO. We then conclude. 



II. THE QUANTUM OPTICAL FREQUENCY COMB OF A SINGLE OPO ABOVE THRESHOLD 



A. Hamiltonian of the system 



We consider the simplest possible case of an OPO with a single, nondegenerate nonlinear interaction. In this case 
the interaction-picture Hamiltonian is 

n 

H^nt = 2ihxP'^{a\a}_^ - a^a^i), (1) 

i=l 

where /3 is the classical (real) and constant (undepleted) pump field (in practice a stable, narrow-linewidth, continuous- 
wave laser) and a±i are the photon annihilation operators of entangled Qmodes ±z, of frequencies 

^±. = Y±(* + ^)a, (2) 

where tjp ~ uji + u^i is the pump frequency (see Fig. 1) and A is the free spectral range of the OPO cavity. 



a. 



a_2 ^-1 Cii a2 flj 



- {2n + 1)A ujp- A 



Ljp + A tjp -I- (2n. + 1)A 



FIG. 1. Optical frequency comb defined by the resonant modes of the OPO cavity, spaced by free spectral range A. The green 
arrow symbolizes the pump field, placed at half its frequency for clarity. 



Below the OPO emission threshold, such a system is known to emit two-mode squeezed fields which demonstrate 
the Einstein-Podolsky- Rosen (EPR) paradox [27, 28]. Above the OPO threshold, the undepleted classical pump 
approximation can still be taken to hold and the generation of EPR states has also been shown to be possible, 
theoretically [ ] and experimentally [;-)0-;-)3] . ['W] 

However, the undepleted classical pump approximation breaks down if the external pump power is increased sig- 
nificantly above threshold. In that case, the pump must be treated as a quantum field in the three-wave mixing 
interaction 

n 

Hint = 'i.ihx ^(paJaL, - p^aia^i), (3) 
i=i 

where p is the annihilation operator of the pump field. Recently, it was predicted [• ! 'i] and experimentally demonstrated 
[30] that the pump field participates in three-way entanglement in this case. Another interesting theoretical analysis 
showed that the signal fields from two OPOs pumped by the same field could become entangled [ : ]. Here, we extend 
this analysis to the QOFC of a single OPO, in which a vast number of different Qmode pairs are already known to be 
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entangled by their parametric downconversion from the pump field [12, 38]. In the undepleted pump approximation, 
all EPR pumps are independent. However, when one considers the OPO well above threshold, there is but a single 
quantum pump field, whose strong (ideally total) depletion entail strong correlations between the EPR fields, since a 
pump photon downconverting into one Qmode pair will necessarily not be downconverted into any other pair. This 
paper posits that this situation should yield multipartite, rather than bipartite, quantum correlations and our goal is 
to ascertain whether they result in multipartite entanglement, which they do. 



III. THE QUANTUM OPTICAL FREQUENCY COMB OF A SINGLE OPO ABOVE THRESHOLD 

As mentioned before, we consider the simplest possible case of an OPO cavity with a single pump mode and a 
single nondegenerate interaction in its nonlinear crystal. Such a crystal implements the Hamiltonian of Eq. (3). We 
assume a single two-mirror standing wave cavity, with one mirror of reflectivity R^^ = 1 for all modes and the other 
an output coupler of R±i = 1 — T±i < 1. Taking into account the vacuum modes Af^ that enter the cavity through its 
output coupler, the input-output theory [39, 40] can be used to derive the equations of motion for the internal cavity 
modes 

di = 2xpali - ki a, + y/2ki Aj" (4) 
a_, = 2xpal - fc_, a_, + y/2k~ A^^^ (5) 



P 



= -2xy^aiQ_j - kpp+ y2fc^ p^n- (6) 



Here k±i = T±i/2r are the loss rates of the cavity mirror for mode i, r being the cavity round trip time. In order to 
solve Eqs. (4-6) we first rewrite field operators as centered fluctuations about their expectation value 

a,; — ai + Sat (7) 
AT = SA^ (8) 
p — Tu + Sp (9) 

Pin = ^in + Spin- (10) 



A. Classical steady-state solutions 

The semiclassical, or mean value, equations follow directly from Eqs. (4-6): 

di = 2xwa*_i - ki ai (11) 
d_i = 2xn7a* - fc-i a^i (12) 

n 

m ~ —2x ^ ^ aiCt^i ~ kpUj + y^2kp Win- (13) 

i 

Now, considering same cavity losses for all signal modes, ki = k-i — ka, the stationary solutions of the two coupled 
equations Eqs. (11-12) for semi-classical mean values are 

2xwa°^*, ^ ka (14) 

2xtna^* =fcaa°.. (15) 

this results in = |q!° J , jro^j — ka/2x and 0? + 0" ^ — = where and 0° are the respective phases of 
aj.j and n7°. For simplicity we take -cuin real and positive therefore based on Eq. (13), cj/' = and (p^ = —(p'^i- The 
stationary solution for the pump field's mean value can then be written as 

n 
i 

where 
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Clearly, the right-hand side of Eq. (16) must be positive and cr = 1 defines the threshold pump field 

th _ / ^ 

2xV 2 ' 



137, 



(18) 



Hence, a = (txJin/tn'Jj)^ is also the pump to threshold power ratio. The classical signal amplitudes are weakly set by 
Eq. (16). 

B. Stability analysis 

The stability of the steady-state solution can be determined by a linearized analysis for small perturbations: 

ai — + Sai (19) 
Wi = vj^ + SvJi. (20) 

Substituting Eqs. (19-20) into Eqs. (11-13), we get 



6aLi = kaSa* + 2xciLi5w — kaSa^ 



(21) 
(22) 



(5n7 = — 2x ^^(a^(5a_j; -I- a^_i6ai) — kpSw, (23) 

i 

where i = l,...,7i, n being the number of signal and idler pairs considered inside the cavity. Defining 6 A = 
(• 



6ai Sa* 6a-i Sa* 



Svu Suj*)'^ We can rewrite Eqs. (21-23) in block matrix form: 





( ] 
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: \ 


/ : 




Sai 






ka 








ka 


• 2xa« 





6 at 




5a* 











ka 








2xa°* 


Sa* 
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~ka 





• 2xa°, 





Sa^i 


dt 


Sa*_, 






ka 








ka 





2xa% 


5a*_i 




dw 
















kp 













v.'. 










-2xa°* . 





kp / 


\8V3* 



(24) 



We derived the eigenvalues of the matrix in Eq. (24) for n = 1,2,3. In all these cases, the eigenvalue sets have the 
following form: 



{A} = {0, . . . , 0, -2fca, . . , -2fca , Ai, Aa, A3, A4}, 

2n-l 2n-l 



where, posing the pump-signal loss ratio n, = kpjka-, 



Al,2 = -2^a ( K ± [k - 8(Vo^ - 1)] 



A: 



3,4 



K + 2 ± J (k + 2)2 - BnVf^ 



(25) 

(26) 
(27) 



Because of the particular symmetry of the problem — namely the block structure of the matrix in Eq. (24), we argue 
that it is reasonable to postulate that Eq. (25) is the general eigenvalue set, Vn, even though a complete inductive 
proof is formally required. For certain initial conditions all 2(2n -I- 1) eigenvalues of the matrix in Eq. (25) can only be 
zero or negative, which ensures the stability of the stationary solution presented in Eq. (16). Equations (26-27) show 
that, as the number of times above threshold a increases, one can always find negative values for Ai_. .^4 by increasing 
the pump-signal loss ratio k, thereby tending towards the doubly resonant OPO, which is always stable. 
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C. Quantum fluctuations 



Now, we rewrite Eqs. (4-6) for the quantum fluctuations around these classical mean values. Notice alj — a^, 

Sai = 2xiSpa,e"l'^ + zuSal,) - kaSai + y^6A^^ (28) 
Sa-i = 2xiSpa,e~"l'^ + mSaj) - kaSa^, + y/2kaSA^i (29) 

n 

Sp = -2x^(ckie''^'^a_i + aie~'''*"6ai) - kpSp+ y/2k^Spin. (30) 

i 

We introduce the generalized field quadrature operators as Qi — {e'^^'a] + e^^^^ai) and Pi — i{e^^'al — e^^^^a.i). Then 
solve these coupled equations we can use the symmetry of the equations in the exchange of the two signal modes and 
introduce the new variables [41] 



Qi-\- ' Qi Q — i 
Qi— ' Qi Q — i 

p,+ = p, + 

The equations of motion for these quadratures are 

5Q^- - 



5Qp 



2kaSQ,- + V^ka SQt 

n 

2x a,6Q,+ - kpSQp + y/2kp SQ^ 



SP,+ = Axa,5Pp - 2ka6P,+ + ^2ka 5P^ 



SP,- = 



'2ka 5Pl'l_ 



SPp = -2x^ a,5P,+ - kpSPp + SP;'' 



(31) 
(32) 

(33) 
(34) 

(35) 
(36) 

(37) 

(38) 
(39) 

(40) 



As seen from Eq. (36) and Eq. (39), the equations for the antisymmetric modes are decoupled from the pump and 
the solutions are, in the frequency domain [41], 

<50°T {n) = - — ^ 6Q% in) (41) 



sp°f{n) = -1 - 



2ka, + ifl 
2ika 



The frequency-domain equations for the symmetric modes are [41] 

inSQ,+ {n) = Axa,6Qp{n) + ^/2ka8Qf+{n) 

imPi+{n) = Axa,sPp{n) - 2kaSPi+{n) + ^/2ka5P^{n) 

n 

inSQpifl) = -2x^a,<5Q»+(r2) - kpSQp{n) + y^5Qf{n) 



(42) 

(43) 
(44) 

(45) 



insPpin) = -2x^aj(5A+(f^) - kp6Pp{n) + ^/WpSP'^in). 



(46) 



These equations can be easily solved for pump and signal- idler pairs. The output quadratures are finally determined 
using input-output relations: 



5QT 



^pout ^ ^2ka5P± - SP. 

^out 



<5g°"' = ^2kp5Qp - SQ'p^ 



2ka5Q± - SQ'^ 

in 
± 



^pout 



'2kp5Pp - SP'^ 



(47) 
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The solutions of Eqs. (43-46) are 



n(^-ikpn + n2-8x^Eja]) ^ 



2 '5Q;,"(f^) (48) 



+^ ^ ^ 1^ ^2fca + zr! (2fc<, + zl7)[(2fc,+il])(fcp + il]) + 8x2Eja']y ^' 

^_^>^a,^P,,,(f^) 



{2ka+in)[(2ka+in){kp + in)+^^ 

5pi"m) f49) 

{2kq + in){kp + in) + 8x2 '^^^ ^^^^ 



''''^ Ea.-jg^.lf^) (50) 



kpn + i(n^- 8k^ J2j a') j^i 



2ka{kp - in) + ikpn + fl^ _ gx^ Ej 



«2 



P ^ ' 2ka{kp + z^!) + ikpn - f72 + 8fc2 Ej a] ^ ^ ' 



ikpn - + 8p E, ^ 



(51) 



2fca(fcp + ifi) + ih 

Substituting the classical solutions Eq. (16) in Eqs. (48-51) and taking £7 = 0, these equations yield: 

SQ°_^^ (52) 
5P°f — ^ oo (53) 
<5Q°T oo (54) 



kakpyj <J rlo^Kp^O K^Kp^ (J 

1. Two-mode squeezing 

In order to quantitatively study squeezing behavior, we assumed equal classical mean values for all pairs. However, 
we can assume any ratio between the classical mean values, as long as they satisfy E" I'^iP = ^"''''i^ = const, 
Eq. (16). Therefore the variances of squeezed and antisqueezed quadratures are 

(56) 
(57) 
(58) 

(59) 



V{Q-^) 


= (iSQ-if) - 


^ 


V{P-^) 




-^• CO 


V{Q+i) 


= {iSQ+^)') - 


— >■ oo 


V{P+i) 


- {iSP+^)') = 


2(a-l) 
na 
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which yields the classic EPR result at threshold (cr = 1), the generation of n independent entangled {i, —i) pairs. 

However, if the OPO is above threshold (cr > 1), the variance of the phase sum, Eq. (59) will increase from zero 
[29, 4J] and will eventually stop being squeezed. It states the well-known fact that for the OPO operating well above 
the emission threshold, twin pairs are not independent EPR pairs due to the pump statistics, unless kp <^ ka [■'^"]- 
This is precisely the mechanism that we rely upon to create multipartite entanglement in this work. We give two 
preliminary examples before turning to the evaluation of precise multimode entanglement criteria. 

2. Multimode squeezing 

We give two examples of squeezed multimode operators which will be useful in the next section. 
First off, specifically combining phase sum operators of two pairs i and j yields 

a, {Pj + P-j ) - (P, + ) 0. (60) 

Even though phase sum operators for each pair become noisier with increasing input pump, this noise can be canceled 
appropriate linear combinations. 

Another interesting example is that of operator X]i=i(-P» + P-i) ~ xPp. In Fig. 2 we plot its variance as a function 

1.0 
0.8 
0.6 
4 
0.2 



1.2 1.4 1.6 1.8 2.(1 



FIG. 2. Plot of nV{Pi + P-i), in red, and V^EILiC^* + ^-0 ~ 2;Pp], in blue, for n = 3 and x = a. When a increases, the 
variance of Pi + P-i increases from zero and approaches the shot noise level. However, the variance of X]r=i(^* ^-0 ~ ^Pp^ 
not squeezed at threshold, subsequently drops from the shot noise level (of value 1 in this graph) and shows squeezing. In this 
particular graph, the minimum of V'E"^i(Pi + P-i) — xPp] occurs at cr = 1.18 but, in general, the value of a for which the 
blue curve reaches its minimum, as well as the value of the minimum itself, depends on the choice of x. 

of cr, for X ^ a. This graph clearly shows that the assumption of the existence of a correlation between all modes 
and the pump is a sensible one. Having established this, we turn to directly testing the existence of multipartite 
entanglement in our system. 

IV. MULTIPARTITE ENTANGLEMENT IN THE OPO WELL ABOVE THRESHOLD 

A. The van Loock-Furusawa inseparability criterion 

The van Loock-Furusawa (vLF) multipartite entanglement criterion [ ] is the multipartite generalization of the 
Duan [ ]-Simon [ ] criterion, itself the continuous-variable formulation of the Peres [-.-jJ-Horodecki [ ] positive 
partial transpose criterion. A density operator is partially separable if and only if it can be written as the convex sum 

P = ^ViPi.ki fe,„ "S) Pj,fc„+i...,fc„, (61) 

i 

where the mode set (fci,...,fcm) is separable from the mode set {km+i, ■ ■ ■ ^ k„). If we define two "entanglement 
witnesses," quadrature operators with arbitrary real parameter sets {hi}i and {gi}i, 

U = hiQi + h2Q2 + h hnQn (62) 

V = giPi + g2P2 + ■ ■ ■ + gnPn, (63) 
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then the separable density operator of Eq. (61) must verify the vLF inequahty [26] 

Vp{u) + Vp{v) > 2 (I VSfci + ■ • • + hk^9kj + \hk^+,9k,„^i +■■■ + hk^gkj) , (64) 

whose violation implies the existence of entanglement between mode set (fc^,. . . ,km) and mode set (fcs,. . . ,fc„). Opera- 
tors u and V were coined variance-based entanglement witnesses for continuous-variable systems by Hyllus and Eisert 
[46], in reference to the original Qbit expectation- value-based entanglement witnesses [4 j, 47, 48]. 



B. Multipartite entanglement in a single, depleted-pump OPO 



In order to demonstrate multipartite entanglement, we examine the conditions for violation of all possible vLF in- 
equalities, corresponding to all possible respective mode partitions such as Eq. (61), and their associated experimental 
regimes. 



1. Pump-signals partition 

We first consider the separability of the sole pump mode from all signal modes. We define ui and vi as 

r, " 



Ml 



v^^Y.^P,+P_i)-xPj,, (66) 
1=1 

with the real parameter x > 0. Based on Eq. (64), the separability of mode p implies 

Si = V{ui) + V{vi) ^ 2 {\h-n9-n + ■■■ + /i-i.g-1 + higi + ■■■ + Kgn\ + \hp9p\) (67) 
^ 2 (2|/ii5i + • • • + /i„5„| + \hpgp\) (68) 

r, « 9 " 

^2(|-^a.H-|— ^a.(-x)|) (69) 
ai ^-^ X0L\ '-^ 

^ i=l ^ i=l 

o " 

(70) 

Here, and in the following, we make the assumption that all classical amplitudes {a^ji are equal, for the sake of 
simplicity. This doesn't lessen the generality of our treatment but makes numerical evaluations easier. Under this 
assumption, we get 

^ 8n. (71) 

Figure 3 displays the maximum violation of the above inequality versus n and cr, for optimized values of the arbitrary 
weight X. As can be seen, there always exist values of (n, a) for which Si—8n is negative, which proves the inseparability 
of the pump mode from the signal modes. Unsurprisingly, entangling a larger number of pairs (n) requires a higher 
pump to threshold power ratio (cr). However, a large pump power (cr) does degrade the inseparability, which might 
just be due to the increasing depletion of the intracavity pump field. 



2. Partition of one (aj , a_j ) EPR pair 



We now study the inseparability of an entangled pair (oj, a^j) from the rest of the signals and the pump. For such 
a partition, we define 

n 2 ^ 

U2 = Qj++^—Qi+-\ -^aiQp (72) 



V2 = Y.'^-P,+ -P,+ ), (73) 
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FIG. 3. Left, sketch of the mode partition studied. Center, plot of the optimum values of a; = Xopt which give maximum 
violation Si (x) — 8n of the vLF inequality, at a given pump to threshold power ratio a and a given number n of mode pairs 
inside the cavity. Right, plot of the maximum vLF inequality violation 5'i(2;opt) — 8n, versus a and n. We took the particular 
case f2 = and = Oj 



and the vLF inequality is 




Figure 4 shows the violation of this inequality for a broad range of parameters. It also demonstrates the necessity of 
applying larger input pump intensity when considering more pairs inside the cavity in order to generate inseparability 
between all pairs, again unsurprisingly. 




FIG. 4. Left, sketch of the mode partition studied. Center, plot of the optimum values of a; = a;opt which give maximum 
violation S2{x) — 8(n — 1) of the vLF inequality, at a given pump to threshold power ratio a and a given number n of mode 
pairs inside the cavity. Right, plot of the maximum vLF inequality violation S2(a;opt) — 8(n — 1), versus a and n. We took the 
particular case = and a,; = ctj 
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3. Pair set partition 

We next turn to partitions {(oi, 61) . . . {ak,bk)}{{ak+i, fofc+i) ■ • • (a„, &„)}. The operators are 

k n 2 

W3=VQ,+ + V ^Qj+ + VajQp (78) 

^-^ ai xai 

/e n 

^3=EE(^'+-^^:' + )' (79) 

and the vLF inequahty is 

5*3 = Viu^i) + V(y-i) ^ 2 (|2/iigi + • • • + 2hkgk\ + |2/ifc+i5rfc+i + • • • + 2/i„g„ + hj>g^\) (80) 
^8fc(n-fc). (81) 

Assuming 1 ^ fc < n, then it is straightforward to show that 8(n — 1) ^ 8fc(n — fc) ^ 2n^. As a consequence, the vLF 
inequahty is automatically violated when vLF inequality 5*2 is, and does not need to be considered separately. 

4. Single signal mode partition 

Finally, we consider partitions of a single signal mode, or several such, all belonging to different EPR pairs («, —i). 
Because the signal mode ai is highly entangled to the mode a_i — and to all other signals by virtue of the preceding — 
the separability test is simple in the present case. It is straightforward to show that the necessary vLF inequality for 
the partition {ai, . . . ,afe}{a_i, . . . ,a_fe, (ofc+i, a.jfc+i)), . . . , (a„,a_„))} , 




54 Q, - + V — (g. - Q-^) + F V + P_, - — (P, + P_,) 



> 2 {\higi H h hkgkl + l^-iff-i + ■ • • h^^g-k + '^H+igk+i + • • •2/i„5„|) 



21-E- + E-1+21E — +E — +0 



>4(n-fc), (82) 

is always violated in the presence of single EPR pair entanglement. This is because the left hand side term of Eq. (82) 
contains EPR nuUificrs [ ], a.k.a. EPR entanglement witnesses, whose squeezed variances tend toward zero. The 
inseparability of any other form of partitions on modes when modes oi and a_i are placed in different partitions, 
can be examined by inequalities similar to 6*4 and with nonzero boundaries. Such inequalities are always violated. 
Therefore, if EPR entanglement is present (the checking of which is a staple of the experimental calibration of a 
regular two-mode squeezer), then the violation of both Si and S2 is a necessary and sufficient condition to mode 
inseparability for all possible partitions in the optical frequency comb of a single OPO. 

C. Entanglement between pairs without considering the pump 

Here we ask the question of the possibility of multipartite entanglement between twins without considering the 
pump field. For that, we rewrite inequality 52 without the pump quadratures: 

5; = !/(«;) + y(«;)^8(n-l) (83) 

with 



4 = Q3+ + E — (84) 

n 

4 = E-^:'+-^'^+- (85) 
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In Fig. 5, we plot the violation of this inequality. As can be seen by comparing with Fig. 4, Inequality 5*2 requires 




1 4 6 e la 
n 



FIG. 5. Plot of the vLF inequality violation 5*2 — 8(n — 1), versus a, the pump to threshold power ratio, and n, the number of 
mode pairs inside the cavity, for Q = and Qi = q^, Vi, j. 

slightly larger a to be violated, compared to 5*2, for small values of n. It shows we need to pump harder (and get 
closer to total depletion) in order to see pure entanglement between twin pairs. The arguments of subsections and 
IV B 3 and IV B 4 may be reused here to complete the inseparability proof. 

D. Optimal entanglement witnesses 

A valid question is whether the entanglement witnesses that were derived here on the basis of physical considerations, 
namely the pump-depletion-induced correlations between EPR pairs, are in fact the optimal entanglement witnesses 
for the system. In other words, do different observables exist that would lead to even more strongly violated vLF 
inequalities? Answering this question should, in turn, inform on what type of entangled state is really generated here. 
The search for optimal entanglement witnesses was addressed by Hyllus and Eisert, using semidefinite programming 
procedures [4()]. While the scope of the present paper is limited to this successful demonstration of large-scale 
entanglement in a simple OPO, the exact nature of the quantum state generated is clearly an interesting followup 
question, on which light can be shed by seeking the optimal entanglement witnesses and checking whether they are 
different from the ones derived above. 



V. CONCLUSION 

We showed that a single OPO operating well above threshold can generate multipartite entanglement in its quantum 
optical frequency comb. We verified the multipartite nature of the entanglement by evaluating the van Loock-Furusawa 
separability criterion over all possible Qmode partitions. We showed that all of these vLF inequalities can be violated, 
for an arbitrary large number of pairs n, simply by increasing the input pump power higher above threshold. 

While the presence of multipartite entanglement in such a simple system is a remarkable feature, it is important to 
keep in mind that the exact type of entanglement that is produced here (GHZ, W, cluster) is difficult to determine. The 
search for optimal entanglement witnesses for this system is a promising approach to illuminate this question. Note 
that previous work has shown that multipartite cluster-state generation, which was experimentally demonstrated in a 
single OPO below threshold [ - ] , should actually fail above threshold [ ' '] . However, the multipartite entanglement that 
we discovered in the simple OPO above threshold would certainly be useful for quantum communication applications, 
such as quantum secret sharing. 
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